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INTRODUCTION  \ , 

In  what  follows  we  present  a mathematical  model^for  pre- 
diction of  linear  elastodynamic  response  of  metal  matrix  compo- 
site penetrators.  Although  it  is  recognized  that  for  the  appli- 
cations envisaged,  the  anelastic  effects  are  the  dominant  ones, 
a study  of  linear  elastic  behavior  is  of  importance  inasmuch 
as  it  yields  information  about  the  composite  properties  required 


for  simulation,  using,  for  example,  the  HELP  code 
material  strength  is  not  neglected.  ^ 


where 


If  a rod  of  a fiber-reinforced  composite  is  subjected 
to  impact  at  one  of  its  ends,  longitudinal  waves  propagate 
through  the  rod  in  such  a manner  that  on  the  macroscopic  scale 
transverse  normal  stresses  and  axial  shear  stress  are  zero 


at  least  to  the  first  degree  of  approximation.  Thus  the  prob- 
lem is  essentially  one  of  uniaxial  stress  rather  than  uniaxial 

' ' [2  3l 

strain.  Since  most  of  the  approximate  models  ' for  longitudi- 
nal wave  propagation  in  the  direction  of  the  fiber  axis  in  a 
unidirectional  fiber  reinforced  composite  treat  the  case  of 
uniaxial  strain  only,  it  is  worthwhile  to  develop  an  analogous 
model  for  the  case  in  which  the  primary  mode  of  propagation 
induces  a state  of  uniaxial  stress. 


The  geometry  of  the  system  we  have  analyzed  is  shown  in 
Figure  1.  It  consists  of  a singly  circular  cylindrical  fiber 
embedded  in  a concentric  circular  cylinder  of  the  matrix  whose 
outer  surface  is  stress  free.  Although  it  would  be  more 
realistic  to  consider  a rod  reinforced  by  multiple  fibers,  our 
model  leads  to  simpler  analysis  and  is  expected  to  provide 
adequate  estimates  for  the  gross  mechanical  response  of  the 
composite  rod  penetrator. 

We  first  derive  a binary  mixture  theory  for  the  rod, 
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using  the  asymptotic  technique  developed  by  Hegemier.  As 
a result  of  this  analysis,  we  are  also  able  to  obtain  (1)  the 
Young's  modulus  for  the  composite  in  the  direction  of  the  fiber 
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axis  and  (2)  the  Poisson's  ratio  These  results,  together 

with  the  results  of  the  corresponding  uniaxial  strain  analy- 
sis, ‘ ' * yield  an  estimate  of  one  more  composite  property  as 
will  be  shown  in  one  of  the  following  sections.  Thus,  three 
of  the  five  moduli  required  for  transversely  isotropic  composites 
can  be  determined  using  the  expressions  derived  in  this  report. 
Although  these  properties  have  previously  been  obtained  else- 
where, our  approach  is  somewhat  different,  as  will  be 
‘obvious  from  the  exposition. 

FORMULATION 

With  the  geometry  and  coordinate  system  shown  in  Figure 
1,  axisymmetric  motion  of  the  rod  is  described  by  the  follow- 
ing equations: 

A.  Conservation  of  Momentum 


x XX 


X xr 


?(“)  . pt“)52  ;(») 
ftA  t r 


B.  Constitutive  Equations 


X+2u  X 


X X+2m  X 


X x+2u 


'F"r 


up/F 
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C.  Interface  Conditions 


a<l> 

xr 


3»>, 

xr  rr 


,U, 


c<2> 

rr 


at  r = r. 


(5) 


D.  Free  Surface  Conditions 


- 0. 


o'2>  - 0 

rr 


at  r 


(6) 


E. 


Appropriate  Boundary  and  Initial  Data 

In  the  above  set  of  equations,  which  defines  a well  posed 
boundary  value  problem  in  (x,r,t)  plane,  we  have  used  the  super- 
script a to  identify  the  fiber  (a=l)  and  the  matrix  {a=2) . The 
object  of  the  subsequent  analysis  is  to  eliminate  the  r coordi- 
nate under  the  assumption  that  composite  microdimension  is 
much  smaller  than  signal  wave  length  J for  problems  of  interest. 
For  this  purpose  we  introduce  the  definitions; 


(m) 


» Typical  signal  wavelength 

= Composite  Young's  modulus  in  the  axial 
direction 


^(m) 


(x,r,t) 
(X,U)  ' 

p'“>  - 

(o  , eo 
' XX ' xr 


Composite  density 


“ ^(m)/P(m)  ' ^(m)  ' (m) 

* 12^ 

» (x/I,  ?/?2» 


(7) 


(X,U) 


(m) 


'°ee’®rr’ 


(m) 

(a) 


^®xx''^xr'®90'°rr^ 
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It  is  noted  here  that  at  this  stage  of  analysis,  the 
appropriate  mixture  modulus  and  density  are  not 

known;  however,  they  will  be  determined  in  a later  section. 

In  terms  of  the  variables  defined  in  Eq.  (7) , the  initial 
boundary  value  problem  (1-6)  reduces  to  the  following  equations. 

A.  Conservation  of  Momentum 


3 + 

X XX 


i 3 /ra<“^\  = p^“^32u^“^ 

r r \ XT  J ^ t X 


(a) 

69  = 2 (a),2„(a) 


B.  Constitutive  Equations 


(8) 


(9) 


, a 

/ XX 


rr 


(a) 


X 

X+2y 

X 


X 

X 

X+2y 


(a) 


(a) 


3 u 

X X 


'Uj./r 


(10) 


/s  u + g^3  u \ 
^ r X X rj 


(a) 


(11) 


C.  Interface  Conditions 


u 


(1) 


rr 


u^2), 

r X 


nt2)  „(1) 

‘'rr  ' ‘'xr 


= u(2) 


,(2) 

xr 
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D.  Free  Surface  Conditions 


(2) 


rr 


0,  a 


(2) 


xr 


0 at  r 


(13) 


MIXTURE  THEORY 


To  derive  an  approximate  theory  that  includes  the  ef- 
fect of  material  inhomogeneity,  we  first  define  averaged  and 
partial  quantities  by 


f (x,t) 


1 


f 


f (x,r, t)  dA 


(14) 


(a) 


f (x,t)  « n 

f‘““>(x,t). 

(15) 

dependent  variable  f,  where 

= [0,r^], 

(16) 

. nrl, 

A«>  . X - A<» 

(17) 

n(“>  . A^‘^>A 

• 

(18) 

(o) 


Obviously,  the  quantity  n denotes  the  volume  fraction  of 
the  a constituent.  If  the  axial  momentum  equation  (8)  is  now 
averaged  according  to  (14)  and  the  relevant  interface  and 
free  surface  conditions  are  utilized,  we  obtain  the  following 
equation  for  conservation  of  average  axial  momentum; 


3 o'lP)  - o<lP>3^u'l“>  . -I 
X XX  “ t X 

3 - b'2p)3?u<2">  - P 


X XX 


t X 


(19) 

(20) 


where  P is  an  interaction  term  reflecting  the  axial  momentum 
transfer  from  the  matrix  to  the  fiber,  and  is  given  by 
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P = 2r^  (x,r,t) 


Equations  (19)  and  (20) , which  are  basic  to  the  approximate 
analysis  of  wave  propagation  in  a composite  rod,  are  to  be 
complemented  by  appropriate  constitutive  relations  for  partial 


stresses  a 


and  the  interaction  term  P.  Such  relations  are 


most  easily  obtained  through  an  asymptotic  scheme  which  shall 
now  be  pursued. 

Asymptotic  Expansions 

2 

The  system  (8-13)  contains  the  parcuneter  e , which  is 
square  of  the  ratio  of  rod  radius  to  the  signal  wavelength. 
For  problems  of  interest  here,  this  quantity  is  much  smaller 
than  unity.  Hence  we  expand  all  the  dependent  variables  in 
a power  series  in  this  parameter;  thus,  we  write 


:,r,t,e)  = ^ 


^(2n) 


for  any  dependent  variable  h.  If  the  expansion  (22)  is  sub- 
stituted into  (8-13)  , a hiearchy  of  models  is  obtained  by 
retaining  the  terms  corresponding  to  different  values  of  n. 

The  simplest  theory  can  be  derived  from  the  zeroth  order  system, 
the  first  equation  of  which  is  obtained  from  (11).  Thus,  we 
have 

3 u^“^  * 0 (23) 

^ *(0) 

so  that 

"x'Jo)  * “x'(0) 

From  the  zeroth  order  expansion  of  the  radial  momentum  equa- 
tion (9)  and  the  necessary  constitutive  equations,  interface 
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and  free  surface  conditions  we  obtain  the  following; 
In-Plane  Stress  Problem 


a -(a)  . 

^r°rr(0) 


\°00(0)/ 

,(2) 
rr(0) 


°rr(0) 


(a)  _ (a)  \ « 

‘^rr(O)  ‘^00(0)1  “ ° 


(a) 


= x'"' 

3 u^“^ 
^x^x(O) 

1 1 

+ 

^X+2u 

1 

/3  u 

' 

X+2uj 

\v-/ 

0 at 

r = 1 

'^rr(O) ' 

“r(i)  = “ 

(2) 

r(0) 

at  r 

= >^1 

(25) 

(a) 

(26) 


(27) 

(28) 


If  the  boundary  value  problem  (25-28)  is  solved  subject  to  the 
condition  that  be  bounded  at  r = 0,  we  obtain 


u 


(1) 

r(0) 


(29) 


(2)  _ A^2)  (2) 

"r(0)  " ^ ^ 


where 


A<2)  . 


- B<2) 


and 


B 


(a) 


e-1 


®aB®x'*x(0) 


(30) 


(31) 


(32) 


7 


R-3401 


The  elements  of  the  matrix  b . in  Eq.  (32)  can  be  expressed  in 

Ot  p 

terms  of  the  material  properties  by  the  following  equations: 


where 


bj^j^  « (1  + n^^^)  /2D 

= n X y VzD  (33) 
bjj  = X«>  (u<'>  * + x'l>)/2D 

D.-[(x<»  + p<l')  {x'2>  * „<2)  (!*„(!>){ 


* (x<2>  ^ 


p<^>) 


Using  (29-31)  we  can  obtain  constitutive  relations  for  the 
partial  stresses  in  the  following  manner.  First,  from  (10) 
we  have: 


XX  (0) 


■ (X«P)'“’  * X'“>  i (35) 


SO  that  on  using  (14,15)  we  obtain 


= rx+2u\^2p)  3 (2a)  (2)  f (2)  . 

xx(0)  ^x"x(0)  ^ ['^r(O) 


(36a) 


°xx(0)  {^+2y)  3x'^x(0)  * [^l'^r(O)  ^^'^1'^^ 

(36b) 

If  (29-32)  are  now  substituted  into  (36) , we  finally  obtain 


f 
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XX  (0) 


^o6  ^x^x(O) 


where 


with  6^g  denoting  the  Kronecker  delta.  It  can  easily  be  shown 
that  the  matrix  is  symmetric,  i.e. , 

^^12  ^ ^21 

axial  shear  problem  and  the  interaction  term 

To  obtain  the  constitutive  relation  for  the  interaction 
term  P,  we  use  the  zeroth  order  expansion  for  the  axial  momen- 
tum equation,  which  is 


a 4. 

^x'^xx(O) 


(■ 


xr  (0), 


- P ^t"x.(0)' 


To  (40)  we  append  the  necessary  constitutive  equations,  inter- 
face and  free  surface  conditions: 


^xr(O) 

a u^“^ 
®r'^x(2) 

+ a u^“^ 

^ ^x'^r(O) 

(41) 

‘^xr(O) 

= rT<2) 

'^xr(O) 

P 

* 2r 

zri 

at  r = rj^ 

(42) 

«(2) 

‘^xr(O) 

= 0 at 

r = 

1 

(43) 

It  is  noted  here  that  in  (42)  we  have  used  the  definition  (21) 
with  0^“^  replaced  by  (^^^^(O)  * 

We  now  use  (35,41)  in  (40)  to  obtain  differential 
(ot) 

equations  for  ^^^{2)'  which  turn  out  to  be 
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(>a) 


r 3 


= (x,t) 


(44) 


(a) 


where  Q represents  terms  in  (40)  which  are  independent 
of  the  radial  coordinate.  We  can,  however,  relate  to 

the  interaction  term  by  integrating  (44)  over  and  using 

(42),  (43)  and  (29)  through  (31).  The  final  result  of 
this  procedure  is 


,(a) 


= -(-1) 


a P _ 5,,  (“);. 


(45) 


n 


Equation  (44)  can  now  be  solved  in  conjunction  with  (41-43) 

and  (45)  to  yield  u^,L.  and,  consequently,  the  displacement 
(a)  xuj  2 

fields  u correct  to  0(e  ).  Thus,  we  obtain 


(1)  _ (1)  2 . ^ 

'"x  ■ "x(0)  ^ ^ H^(x,t)  + ^ “(1)  117 

\ V n 

') 


- 2 


(46a) 


m 


7^ 


n 


(1) 


n<^^  - £n  r^  + r^ 


- n^^> 


- n^^^  £n  r^  + r^ 


y I (n' 


(n^^J  iln  n^^)  - n^^^ 


- n(^>  £n  r2) 


(46b) 


where  Hj^(x,t)  is  an  undetermined  function.  To  proceed  further, 
(46)  is  averaged  according  to  (14)  to  yield 
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, (la)  _ , (1)  . ^2„  . ^2  n 
u..  = U.  ...  + E H,  + c — 


(1) 


'x(O) 


7^^ 


2 d B 

X 


(47a) 


(2a)  (1)  ^ 2 

U = U + E 

X X (0) 


„ +.2fE(_l 1 /tn 


g pC2)/n^^^  . ^(1)  . 

®x®  I -JIT  ~ — 


) 


+ y 3 
2 X 


(47b) 


From  (47)  we  obtain  the  desired  expression  for  P in  terms  of 
average  displacements,  i.e., 


P * 


u(2a)  . u(la)  2 

X X 


eS 


a=l 


(48) 


where 


Y = Q 


(1)  „(2) 


1 /,.-(!)  , 2 In  n^^’  \ 


(49a) 


(1) 


1_ 

Ay 


n<»  ( 


2 In  n‘‘’\ 


(49b) 


(2)  . 2n^^^  In  n^^^ 
+ ^2^; 


n 


(49c) 


In  order  to  complete  the  formulation  we  replace  i*'  ^32) 


(aa) 


by  the  corresponding  averages,  i.e.,  bV  ; this  procedure 

furnishes 
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(2a)  .(la)  2 


u - u 

X X 


e Y 


(a)  " X 


where 


In  a similar  fashion,  we  drop  the  subscript  zero  from  the 
variables  in  (37),  thus,  obtaining. 


EC  o 9 u' 
aB  X > 


COMPOSITE  PROPERTIES 

The  foregoing  analysis  completes  the  construction  of 
an  approximate  theory  for  longitudinal  wave  propagation  in  a 
composite  rod,  which  is  given  by  the  equations  (19)-(20),  (33)- 
(34)  and  (49)  through  (52).  To  obtain  the  mixture  Young's 
modulus  from  these  equations  we  eliminate  from  them  all  of  the 
variables  except  , to  obtain 

[(=11  * ^<=12  ^ -=22)  ’t 

+ 0 (€^)]  - 0 (53) 


Equation  (53)  suggests  the  definition  of  mixture  density  ac- 
cording to 


r 


. s 
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We  choose  the  mixture  modulus  such  that  the  coefficient  of  9^ 
in  (53)  is  unity.  This  choice,  together  with  (7,  33,  34  and 
38),  furnishes 


(m) 


(x'l'  . . (x«'  . 2m»i)  „'2)  . eVd 


(54) 


where 


E*  - 


x'i>  lx»i  (l.„ 


)| 

+ X‘2)„(2)„U)]  ^ „(2)y(2)  |^(1)-(2)„(1) 

. x<2>  (x<i>  . * ?«)„«>)] 

- (x<l>  . [x«'  . 7(2)  (i  * „(!)) 

* 7<2),>(2)  (r(2)  . 7(2)) 


(55) 


To  obtain  one  of  the  mixture  Poisson's  ratios,  we 
I first  note  that  the  stress-strain  relations  for  a transversely 

i isotropic  material  are  given,  in  part,  by 

( 


, S. 
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On  the  other  hand,  for  the  composite  cylinder,  we  have,  from 
f39)-(31),  that 


0=1  6=1 


(58) 


If  we  now  set 


a 11  « a 11^2)  a ,, 

0 XI  /A\  ® B XI  /r\\  **  B Xi 
X X (0)  X X (0)  X X 


(59) 


we  obtain 


'^rlo) 


t t 

0=1  6=1 


3u 


3x 


(60) 


A comparison  of  (57)  and  (60)  suggests  the  definition  of  equi- 


valent Poisson's  ratio  for  the  composite  by 


2 2 

V =-  y 

xr 

0=1  6=1 


'o6 


(61) 


Finally,  on  using  (33-34)  and  reverting  to  dimensional  material 
properties,  we  obtain 


xr 


(x'2'  + 2J<2))  „<1>  + x«'  (x'l'  * ?<1> 

;:«>)  n«']/2D 


(62) 


where  D is  given  by  (55). 

As  was  mentioned  in  the  introduction,  we  can  obtain  one 


more  combination  of  material  properties  by  using  the  results 
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of  the  corresponding  uniaxial  strain  analysis.  To  do 
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(n\ 

where  v'  denotes  the  Poisson's  ratio  for  the  a constituent, 
we  can  use  (67)  to  obtain  the  transverse  Young's  modulus,  E^, 
of  the  composite.  Thus  the  only  composite  property  that  re- 
mains to  be  determined  is  the  axial  shear  modulus  to  be  used 
in 


a 


rx 


2u 

X 


e 


rx 


(69) 


We  can  calculate  from  the  expression  given  in  [6],  i.e., 
from 


:-(2)  ir<l>  (!.„»>).  u«'n«' 

“ Fav.(2)  . 


+ U 


d+n'*') 


(70) 


Thus,  all  the  elastic  moduli  of  the  ccxnposite  can  be  calculated 
from  the  expressions  (54,  62,  67,  68  and  70). 


NUMERICAL  RESULTS 

Using  the  expressions  given  here  for  the  composite 
properties,  we  have  calculated  the  elastic  moduli  and  Poisson's 
ratios  for  five  sets  of  material  combinations  as  a function  of 
the  fiber  volume  fraction.  The  constituent  properties  used 
for  these  calculations  are  shown  in  Table  1,  The  calculated 
composite  properties  are  given  in  Tables  2 through  6.  From 
these  tables  it  can  easily  be  concluded  that  calculation  of  the 
Poisson's  ratio  by  using  the  rule  of  mixtures,  which  was 
used  for  computing  v^,  does  not  entail  any  significant  loss 
of  accuracy. 

From  the  results  shown  in  these  tables,  various  elas- 
tic wave  speeds  can  very  easily  be  calculated  since  the 
density  of  the  composite  is  obtained  simply  by  the  rule  of 
mixtures.  Although  the  dispersion  curve  for  time  harmonic 
longitudinal  waves  in  the  composite  rod  can  also  be  computed 
using  the  theory  derived  in  this  report,  we  have  chosen  not 
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to  do  so  since  dispersive  phenomena  are  not  expected  to  be 
important  factors  in  materials  selection  for  metal  matrix 
composite  penetrators. 

CONCLUDING  REMARKS 

A binary  mixture  theory  has  been  derived  for  propagation 
of  longitudinal  waves  in  a composite  rod,  which  is  proposed 
a simple  model  for  metal  matrix  composite  penetrators.  Theory 
includes  the  effect  of  material  inhomogeneity  and  exhibits 
dispersive  character  typical  of  composite  materials. 

Expressions  for  computing  the  composite  material  proper- 
ties have  been  presented,  and  have  been  used  to  calculate  the 
composite  moduli  for  five  sets  of  material  combinations  as  a 
function  of  fiber  volxime  fraction. 
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TABLE  1 

CONSTITUENT  ELASTIC  PROPERTIES 


Young ' s 
Modulus 


Poisson's 


(GPa) 

Ratio 

Aluminum  6061-T6 

69 

0.33 

Steel  1010 

200 

0.29 

Lead 

14 

0.43 

Sic  on  W Substrate 

427 

0.19 

Aluminum  FP 

345 

0.27 

Thorne 1 50  Gr 

393 

0.20 

Tungsten  Wire 

407 

0.28 

R-3401 


TABLE  2 

COMPOSITE  ELASTIC  PROPERTIES 
ALUMINA  FP  FIBERS  IN  ALUMINUM  MATRIX 


( 


I 


Fiber 

Volume 

Fraction 

Young's  Moduli 
E E 

*(GPa)  ^ 

0.2 

124 

97 

0.4 

180 

129 

0.6 

235 

171 

O 

• 

CO 

290 

235 

Axial 

Shear 


Poisson' 

'^r 

s Ratio 

''xr 

Modulus 

Wx 

(GPa) 

0.318 

0.316 

34 

0.306 

0.303 

45 

0.294 

0.291 

62 

0.282 

0.280 

88 

► 


21 


n 


1 


f 


I 

r 


1 

i 

! 


t 

z.  • 

f •* 

S '*• 

t 


i 


i 

t 

! 


i 


i 
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TABLE  3 


COMPOSITE  ELASTIC  PROPERTIES 
THORNEL  50  Gr  FIBERS  IN  ALUMINUM  MATRIX 


Fiber 

Volume 

Fraction 

Young's  Moduli 

E E 

^(GPa)  ^ 

Poisson ' 

''r 

s Ratio 

V 

xr 

Axial 

Shear 

Modulus 

Vx 

(GPa) 

o 

• 

to 

134 

103 

0.304 

0.300 

""  35 

0.4 

199 

140 

0.278 

0.272 

47 

0.6 

264 

189 

0.252 

0.246 

66 

o 

• 

CD 

328 

264 

0.226 

0.222 

98 

1 


2 


I 

1 

t 

I 

I 

i 

i 


I 

1 


r 


t 


22 
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TABLE  4 

COMPOSITE  ELASTIC  PROPERTIES 
Sic  (ON  W SUBSTRATE)  FIBERS  IN  STEEL  MATRIX 


Fiber 

Volume 

Fraction 

Young's  Moduli 

E E 

’'(GPa)  ^ 

Poisson 
■ '^r 

's  Ratio 

V 

xr 

Axial 

Shear 

Modulus 

Vx 

(GPa) 

CM 

• 

O 

246 

237 

0.270 

0.268 

91 

0.4 

291 

277 

0.250 

0.247 

107 

0.6 

336 

321 

0.230 

0.227 

126 

o 

• 

00 

382 

370 

0.210 

0.208 

150 

I 
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I 


I 
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TABLE  5 

COMPOSITE  ELASTIC  PROPERTIES 
TUNGSTEN  WIRE  FIBERS  IN  STEEL  MATRIX 


Fiber  Young's  Moduli 


Volume 

Fraction 

E E 

*(GPa)  ^ 

CM 

• 

O 

241 

230 

0.4 

283 

2S4 

0.6 

324 

303 

o 

• 

00 

366 

350 

Poisson ' 

''r 

's  Ratio 

V 

xr 

Axial 

Shear 

Modulus 

Wx 

(GPa) 

0.288 

0.288 

89 

0.286 

0.286 

102 

0.284 

0.284 

118 

0.282 

0.282 

136 

24 
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TABLE  6 


COMPOSITE  ELASTIC  PROPERTIES 


ALUMINA  FP  FIBERS  IN  LEAD  MATRIX 


i 


Fiber 

Volume 

Fraction 

Young's  Moduli 

E E 

^(GPa)  ^ 

Poisson' 

''r 

s Ratio 

V 

xr 

Axial 

Shear 

Modulus 

Vx 

(GPa) 

0.2 

80 

39 

0.398 

0.395 

7 

0.4 

146 

59 

0.366 

0.362 

11 

0.6 

213 

90 

0.334 

0.330 

17 

0.8 

278 

150 

0.302 

0.299 

33 

{ 

, -v 


